NASA  Contractor  Report  189607 
ICASE  Report  No.  92-26 


ICASE 


j 


CROSSFLOW  EFFECTS  ON  THE  GROWTH  RATE 
OF  INVISCID  GORTLER  VORTICES 
IN  A  HYPERSONIC  BOUNDARY  LAYER 


Yibin  Fu 
Philip  Hall 


DTIC 

ELECTE 
AUG  13  1992 


A 


Contract  No.  NAS  1-1 8605 
June  1992 


Institute  for  Computer  Applications  in  Science  and  Engineering 
NASA  Langley  Research  Center 
Hampton,  Virginia  23665-5225 

Operated  by  the  Universities  Space  Research  Association 


NASA 

National  Aeronautics  and 
Space  Administration 

Langley  Research  Center 

Hampton,  Virginia  23665-5225 


9S  8  11  083 


iiiiiiii 


This  document  has  been  approved 
for  public  release  and  sale;  its 
distribution  is  unlimited. 


CROSSFLOW  EFFECTS  ON  THE  GROWTH  RATE  OF  INVISCID 


GORTLER  VORTICES  IN  A  HYPERSONIC 
BOUNDARY  LAYER 


Yibin  Fu  and  Philip  Hall* 

Department  of  Mathematics 
University  of  Manchester 
Manchester  M13  9PL 
UNITED  KINGDOM 

ABSTRACT 

The  effects  of  crossflow  on  the  growth  rate  of  inviscid  Gortler  vortices  in  a  hypersonic 
boundary  layer  with  pressure  gradient  are  studied  in  this  paper.  Attention  is  focussed  on  the 
inviscid  mode  trapped  in  the  temperature  adjustment  layer;  this  mode  has  greater  growth 
rate  than  any  other  mode.  The  eigenvalue  problem  which  governs  the  relationship  between 
the  growth  rate,  the  crossflow  amplitude  and  the  wavenumber  is  solved  numerically,  and  the 
results  are  then  used  to  clarify  the  effects  of  crossflow  on  the  growth  rate  of  inviscid  Gortler 
vortices.  It  is  shown  that  crossflow  effects  on  Gortler  vortices  are  fundamentally  different  for 
incompressible  and  hypersonic  flows.  The  neutral  mode  eigenvalue  problem  is  found  to  have 
an  exact  solution,  and  as  a  by-product,  we  have  also  found  the  exact  solution  to  a  neutral 
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1  Introduction 


This  paper  continues  our  studies  on  the  Gortler  instability  mechanism  in  hypersonic  bound¬ 
ary  layers.  For  a  discussion  of  the  motivation  for  such  studies  and  background  information, 
the  reader  is  referred  to  our  previous  papers,  Hall  and  Fu  (1989)  and  Fu,  Hall  and  Blackaby 
(1992).  Both  are  concerned  with  linear  stability  properties  of  Gortler  vortices,  but  the  first 
studies  a  Chapman’s  law  fluid  and  the  second  is  devoted  to  a  Sutherland’s  law  fluid.  A  com¬ 
mon  feature  between  the  results  for  these  two  types  of  fluids  is  that  the  most  rapidly  growing 
Gortler  vortices  in  a  hypersonic  boundary  layer  are  those  trapped  in  the  temperature  adjust¬ 
ment  layer  sitting  at  the  edge  of  the  boundary  layer  over  which  the  basic  state  temperature 
decays  from  its  0(M 2)  values  near  the  wall  to  its  free  stream  value,  here  M  is  the  free  stream 
Mach  number.  The  major  differences  when  these  two  different  viscosity  laws  are  used  are  as 
follows.  Firstly,  using  Chapman’s  viscosity  law  predicts  an  0(M 2)  Gortler  number,  whilst 
using  Sutherland’s  viscosity  law  predicts  a  smaller  0(A/3/2)  Gortler  number.  Secondly,  the 
evolution  of  Gortler  vortices  in  Chapman  fluids  is  less  strongly  affected  by  nonparallel  effects 
than  those  in  Sutherland  fluids.  This  is  due  to  the  fact  that  when  Chapman’s  law  is  used,  the 
temperature  adjustment  layer  is  logarithmically  thin  whilst  when  Sutherland’s  law  is  used, 
the  adjustment  layer  is  much  thicker.  Certainly,  for  a  hypersonic  boundary  layer  across 
which  the  temperature  varies  significantly,  it  is  more  appropriate  to  use  Sutherland’s  law 
to  model  the  viscosity  of  the  fluids.  Thus,  in  our  later  papers,  Fu  and  Hall  (1992a,  1992b) 
which  study  nonlinear  development  and  secondary  instabilities,  Sutherland’s  viscosity  law 
has  been  consistently  used.  In  Fu  and  Hall  (1992a),  we  have  shown  how  Gortler  vortices  may 
grow  in  the  neighbourhood  of  the  neutrally  stable  position  to  form  a  large  amplitude  vortex 
structure  which  consists  of  a  core  region  bounded  by  two  viscous  transition  layers.  Once 
such  a  structure  has  been  established,  various  secondary  instabilities  in  the  form  of  traveling 
waves  may  be  triggered  off.  There  are  at  least  two  such  secondary  instabilities  which  have 
been  observed  in  experiments.  The  first  is  called  the  wavy  vortex  mode  of  instability  which 
makes  the  vortex  boundaries  wavy  and  which  is  located  in  the  two  viscous  transition  layers. 
This  secondary  instability  has  been  studied  in  Fu  and  Hall  (1992a).  The  other  type  is  an 
inviscid  instability  which  can  exist  because  of  the  inflexional  nature  of  the  steady  state.  In 
Fu  and  Hall  (1992b),  we  have  shown  how  the  change  in  the  basic  flow  caused  by  the  large 
mean  flow  correction  of  Gortler  vortices  would  affect  the  growth  rate  of  Rayleigh  instability. 

All  of  our  previous  studies  in  the  hypersonic  context  have  been  concerned  with  problems 
which  arise  when  Gortler  vortices  occur  in  two-dimensional  boundary  layers,  but  in  many 
practical  situations,  the  basic  boundary  layer  may  be  three-dimensional.  Boundary  layer 
flows  over  a  three-dimensional  obstacle  or  over  a  turbine  blade  are  such  examples.  A  more 
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significant  example  is  the  boundary  layer  flow  over  the  two  areas  of  concave  curvature  on  the 
lower  side  of  a  modern  laminar  flow  airfoil  which  is  fully  three-dimensional  when  the  wing 
is  swept.  Obviously,  in  order  to  effectively  control  or  delay  transition,  the  role  played  by  a 
crossflow  component  in  the  basic  flow  should  be  fully  understood. 

This  paper  may  be  considered  as  extending  the  Hall  (1985)  and  Bassom  and  Hall  (1991) 
studies  on  crossflow  effects  in  incompressible  flows  to  hypersonic  flows.  It  is  also  related  to 
the  recent  work  by  Dando  (1992)  on  crossflow  effects  in  compressible  flows.  All  of  these 
studies  have  shown  that  crossflow  has  a  stabilizing  effect  on  Gortler  vortices,  and  more 
importantly,  Gortler  vortices  may  be  completely  stabilized  by  the  three-dimensional  nature 
of  the  basic  state.  One  objective  of  the  present  paper  is  to  give  a  quantitative  description  of 
the  stabilizing  effect  of  a  crossflow  on  Gortler  vortices  in  hypersonic  boundary  layers. 

As  has  been  observed  before  by  Hall  (1985)  for  incompressible  flows,  the  case  when  the 
crossflow  velocity  component  is  a  constant  multiplier  of  the  streamwise  velocity  component 
is  degenerate  in  the  sense  that  by  a  suitable  transformation,  the  problem  can  be  reduced  to 
a  Gortler  problem  corresponding  to  a  two-dimensional  basic  state.  This  degeneracy  happens 
when  there  is  no  pressure  gradient  in  the  inviscid  flow.  Thus  in  order  to  reveal  the  genuine 
implications  of  having  a  crossflow,  we  avoid  this  degeneracy  by  assuming  in  this  paper 
that  there  exists  a  pressure  gradient  in  the  inviscid  flow.  However,  as  has  been  shown 
by  Stewartson  (1964),  when  there  is  a  pressure  gradient  in  the  boundary  layer,  analytical 
solutions  to  the  boundary  layer  equations  are  possible  only  if  we  make  three  assumptions, 
one  of  which  is  that  the  viscosity  should  be  given  by  Chapman’s  law.  In  order  to  make  the 
necessary  analytical  progress,  we  shall  assume  that  this  is  the  case.  Thus  the  present  paper 
is  more  relevant  to  Hall  and  Fu  (1989). 

An  important  property  of  a  hypersonic  boundary  layer  is  that  the  chordwise  acceleration 
of  a  fixed  fluid  particle  in  the  basic  state  is  asymptotically  large  and  in  the  opposite  direction 
to  the  centrifugal  acceleration  due  to  the  curvature  of  the  wall.  This  negative  chordwise 
acceleration  has  been  interpreted  in  our  previous  papers  as  effectively  producing  a  negative 
curvature  (or  equivalently  a  negative  Gortler  number).  Thus  the  curvature  of  the  wall  should 
be  large  enough  to  overcome  this  negative  curvature  due  to  the  basic  state  if  Gortler  vortices 
are  to  exist  at  all.  In  this  paper,  we  assume  that  the  Gortler  number  is  of  the  same  order 
as  that  of  the  effective  negative  Gortler  number  due  to  the  basic  state  curvature.  Then  the 
large  resultant  Gortler  number  drives  a  large  growth  rate  and  as  a  result,  the  eigenvalue 
problem  determining  the  growth  rate  is  essentially  inviscid.  The  main  aim  of  this  paper  is 
to  determine  how  the  existence  of  a  crossflow  would  modify  the  growth  rate  of  these  inviscid 
Gortler  vortices. 

The  rest  of  this  paper  is  organized  as  follows.  As  a  preparation  to  later  analysis,  we  first 
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discuss  in  the  next  section  the  solutions  to  the  three-dimensional  basic  state  when  there  is 
a  pressure  gradient.  In  section  3  we  formulate  the  general  linear  Gortler  instability  problem 
which  is  then  specialized  in  section  4  to  the  inviscid  problem  mentioned  above.  Sections 
4  and  5  are  devoted  to  the  non-neutral  and  neutral  solutions  respectively,  and  in  the  final 
section  we  discuss  our  results  and  make  some  conclusions. 


2  The  basic  state 


Consider  a  hypersonic  boundary  layer  over  a  rigid  wall  of  variable  curvature  (\ / A)k{x* / L), 
where  L  is  a  typical  streamwise  length  scale  and  A  is  a  lengthscale  characterizing  the  radius 
of  curvature  of  the  wall.  We  choose  a  curvilinear  coordinate  system  with  x * 

measuring  distance  along  the  wall,  y *  perpendicular  to  the  wall  and  z*  in  the  spanwise 
direction.  To  fix  ideas,  we  may  assume  that  the  rigid  wall  is  the  surface  of  an  infinite 
cylinder;  the  z-axis  is  then  aligned  with  a  generator  of  the  cylinder  (see  Fig.l).  The  velocity 
components  are  denoted  by  (u*,  v*,  w*)  and  density,  temperature  and  viscosity  by  />*,  T*  and 
p.*  respectively.  The  free  stream  values  of  these  quantities  will  be  signified  by  a  subscript  oo. 
The  Reynolds  number  R,  the  curvature  parameter  A  and  the  Gortler  number  G  are  defined 
by 

A  =  4,  G  =  2R1/2A.  (2.1) 

Moo  A 

We  assume  that  the  Reynolds  number  is  large,  while  A  is  taken  to  be  small.  In  the  context 
of  incompressible  flows,  the  characteristic  values  of  G  for  instability  are  0(1);  whilst  in 
hypersonic  flows  G  is  0(M2)  when  Chapman’s  viscosity  law  is  used  (see  Hall  and  Fu  (1989)), 
and  is  0(M3^2)  when  Sutherland’s  viscosity  law  is  used  (see  Fu,  Hall  and  Blackaby  (1992)).  In 
the  following  analysis,  coordinates  (a:*,  y *,  z*)  are  scaled  on  (L,  R~'/2L,  R~l^2L),  the  velocity 
(um,v*,w*)  is  scaled  on  (u^,, /2_1/,2u^0,  R~1/2u^0)  and  other  quantities  such  as  p*,T *,  and  p* 
are  scaled  on  their  free  stream  values  with  the  only  exception  that  the  pressure  p*  is  scaled 
on  and  the  bulk  viscosity  A*  is  scaled  on  p^.  All  dimensionless  quantities  will  be 

denoted  by  the  same  letters  without  a  superscript  *.  For  an  ideal  gas  without  dissociation, 
the  continuity,  Navier-Stokes,  energy  equations  and  the  equation  of  state  are,  to  leading 
order,  given  by 

t  +  a^>  =  0'  (2'2) 

Du  dp  d  ,  du.  d  .  du 
9 ~Dt  ~  ~di  +  +  ITz^lh*'  (2'3) 

Dv  1  2  dp  d  2  d  dvp 

p(~Di  +  2GKU^~~Rty  +  di{(  ~  Z^dTp)  +  d70^^ 


(pvp)  =  0, 

(2.2) 

du .  d  du 

^dy’+  dz^dz’' 

(2.3) 
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(2.4) 


d  .  dv.  d  .  dv 

Dw  dp  dj  2  dvp\  d  dvp 

p  Dt  -  fla2  +  d2\(A  3^8xsJ  +  3x/'‘ ' 

+  5(^)  +  S°,S)’  (2-5) 

,1  d  ,  dT.  1  d  ,  dT,  /ft  _x 

+  <rdt//idJ^  +  crdz^dz^  ^ 

7M2p  =  pT.  (2.7) 

• 

Here  terms  of  relative  order  R~l  have  been  neglected  and  we  have  used  a  mixed  notation  in 
which  (i?i,  V2,  V3)  is  identified  with  (u,u,  to)  and  (xi,£2>£3)  with  ( x,y,z ).  Repeated  suffices 
P  signify  summation  from  1  to  3.  The  constants  7,  M  and  a  are  in  turn  the  ratio  of  specific 
heats,  the  free  stream  Mach  number  and  the  Prandtl  number  defined  by 


~  =  m2  —  -u_°°  - 

7  <*_’  <’ 


PooCpa 


where  S  is  the  gas  constant,  k  is  the  coefficient  of  heat  conduction,  and  a 00  =  is 

the  speed  of  sound  in  the  free  stream.  In  equations  (2.3)— (2.6),  the  operator  D/Dt  is  the 
material  derivative  and  is  given  by 


D  d  6  d 

Dt  ox  ay  dz 


The  basic  state  is  given  by 


(u,  v,  w)  =  (u(x,y),  v(x,y),  w(x,y)),  T  =  T{x,y), 

P  =  p{x,y),  p  =  p{x,y),  (2.8) 

where  we  have  allowed  for  a  non-zero  crossflow  component  w(x,  y).  Such  a  three-dimensional 
boundary  layer  may  easily  be  obtained  if,  for  example,  the  angle  of  yaw  of  the  cylinder  A 
shown  in  Fig.l  is  non-zero. 

By  substituting  (2.8)  into  the  governing  equations  (2.2)— (2.7),  we  find  that  the  basic 
state  satisfies 

= °-  <2-9> 


_  ( _du  du\  dP\  d  ( _du\ 

P\*T*+VYy)=-i;  +  Ty\%)' 


(2.10) 


4 


(2.11) 


f _dw  _dw\  d  ( _dw\ 

\  dx  +V  dy)  dy  \  dy)  ' 


2dPx_ 

dxU 

( 

W  ’ 

(2.12) 

(2.13) 

1  d  (  df  \  . 

7  M2PX  =  pT, 

where  P\  is  the  inviscid  solution  for  the  pressure  evaluated  at  y  =  0.  The  equation  obtained 
from  the  y-momentum  equation  (2.4)  is  simply  dp/dy  =  0  which  implies  that  p  =  p(x)  =  Px. 
From  the  inviscid  theory,  we  have  the  relation 


dPx  „  dU , 


(2.14) 


where  p\  and  U\  have  the  same  meaning  as  Px.  If  we  introduce  7\,  the  inviscid  solution  for 
the  temperature  evaluated  at  y  =  0,  then  equation  (2.13)  can  also  be  written  as 


pT  =  piTi. 


(2.15) 


In  this  paper,  we  assume  that  dP^/dx  /  0,  i.e.  there  exists  a  pressure  gradient  in  the 
inviscid  flow.  If  dPi/di  =  0,  then  (2.10)  and  (2.11)  together  with  the  boundary  conditions 
imply  that  w  is  a  constant  multiplier  of  u,  in  which  case  the  original  three-dimensional  basic 
state  can  be  made  two-dimensional  by  rotating  the  coordinate  system  about  the  y-axis  by 
an  angle  given  by  tan_1(ii>/u),  see  Hall  (1985). 

We  now  turn  to  the  solution  of  the  basic  state  equations.  The  reader  is  referred  to  the 
book  by  Stewartson  (1964)  for  a  detailed  discussion  of  these  equations.  Here  we  shall  just 
give  an  outline  of  the  necessary  procedures  and  the  main  results. 

The  immediate  consequence  of  including  a  pressure  gradient  is  that  similarity  solutions 
are  now  only  possible  if  we  use  Chapman’s  viscosity  law 


p  =  C{x)T, 


(2.16) 


assume  that  a  —  1  and  make  the  homenergic  assumption.  (We  recall  that  analytical  similar¬ 
ity  solutions  are  possible  without  making  any  of  the  above  assumptions  if  there  is  no  pressure 
gradient).  In  this  paper,  we  shall  adopt  these  assumptions  to  make  the  necessary  analytical 
progress.  Then  the  first  step  in  the  solution  is  to  determine  the  temperature  T  explicitly.  To 
this  end,  we  first  eliminate  dP\/dx  from  (2.10)  and  (2.12)  to  obtain  the  following  equation: 


(2.17) 
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where 


I  —  T  +  -(7  -  1)A/2u2. 


(2.18) 


On  comparing  (2.17)  with  (2.11),  we  deduce  that  the  general  solution  to  (2.17)  can  be  written 
as 

/  =  —Aw  -f  const., 

where  A  is  an  arbitrary  constant  to  be  fixed  by  boundary  conditions.  Hence 

T  +  ^(7  —  1  )M2u2  4-  Aw  =  const. 


=  T,  +  -(1  -  \)M7Uf  +  AWU 


(2.19) 


where  Wx  is  the  inviscid  crossflow  velocity  component  evaluated  at  y  =  0  and  is  given  by 

...  Uoo  Sin  A  9  fry 

=  TFvhF  =  R tanA’ 

oo 

since  u ^  =  Uoo  cos  A  (see  Fig.  1 )  and  W\  has  been  scaled  on  R~xl2u*00.  Guided  by  previous 
results  for  incompressible  flows,  we  assume  that  the  yaw  angle  A  is  of  order  R~1^2  (this  is 
the  order  at  which  crossflow  begins  to  affect  the  growth  of  Gortler  vortices),  so  that  W\  is 
an  order  one  constant  under  the  large  Reynolds  number  limit.  We  see  that  the  the  size  of 
W\  can  be  adjusted  by  varying  A. 

By  the  homenergic  assumption,  the  total  energy  (scaled  by  CpToo)  in  the  inviscid  flow  is 
a  constant  and  so 

T,  +  i(7  -  1  )M1u;  =  1  +  i(7  -  1  )M\  (2.20) 

where  the  left  and  right  hand  sides  are  the  scaled  energy  evaluated  at  y  =  0  and  y  =  oo, 
respectively.  With  the  use  of  (2.20),  (2.19)  can  be  written  as 


f  +  ^( 7  -  1)A/V  4  A{w  -  Wx)  =  1  4  ^(7  -  1  )M2. 


(2.21) 


If  we  denote  the  temperature  at  the  wall  by  Tw,  then  (2.21)  together  with  u(0)  =  w(0)  =  0 
implies  that 

*  =  (2.22, 

and  so  (2.21)  becomes 


T  =  1  +  i(7  -  1)M!(!  -  uJ)  +  (1  -  -  1  -  i(7  -  1)MJ}  . 


(2.23) 


In  the  special  case  when  the  wall  is  insulated  ( dT/dy  =  0  at  y  =  0),  it  is  easy  to  show  with 
the  aid  of  (2.23)  that  the  wall  temperature  must  necessarily  be  given  by 


■^“'Iwall  insulated  —  ^  +  o(')'  0^  — 


(2.24) 
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where  Tr  is  usually  referred  to  as  the  recovery  temperature.  Once  we  know  the  expression 
for  T,  we  can  now  solve  (2.9)— (2.13)  by  introducing  a  variable  Y  by 


=  J  Pdy 


(2.25) 


and  expressing  u  and  v  in  terms  of  a  function  tp(x,  Y)  as 

fi  =  W'  5  =  ~Ytx  +  (226) 

Then  the  continuity  equation  (2.9)  is  satisfied  and  (2.10),  (2.11)  become 

dYdxdY  dxdY*  dx  +  C{xY<M  P' gyi’  <2'27) 

dtp  dw  dtp  dw  ,n  d2w 

W~dx  ~  ~dx~dY  ~  P'dY*’  ^2'28^ 

Before  we  could  look  for  similarity  solutions  for  the  above  equations,  we  have  to  eliminate 
the  appearance  of  u  through  T  on  the  right  hand  side  of  (2.27).  This  can  be  achieved  by 
introducing  two  new  variable  ( X ,  Y)  through 


(2.27) 


(2.28) 


X  =  fiM2Pt(x)Jr^T)C(x)dx,  Y  =  JUx)Y. 


In  terms  of  X  and  Y ,  equations  (2.27)  and  (2.28)  can  be  written  as 


dtp  d2tp  dtp  d2tp  d3tp  dV\ 
dYdXdY  ~  dXdY *  =  dY*  +  V'~dX 


dtp  dw  dtp  dw  d2w 
dYdX  ~  dXdY  =  dY2' 
after  (2.15),  (2.13),  (2.14)  and  (2.20)  have  been  used.  Here 

v'w  “  & 


(2.29) 


(2.30) 


(2.31) 


(2.32) 


Equations  (2.30)  and  (2.31)  are  to  be  solved  subject  to  the  following  boundary  conditions: 


tp  =  —  =  0,  w  =  0  at  K  =  0, 
dy 


Vi,  w  —*W\  as  F  — ►  oo. 


(2.33) 


The  above  boundary  value  problem  admits  similarity  solutions  of  the  form 


tp  =  a{X)f(r}),  w  =  Wigfr),  r)  = 


b(xy 


(2.34) 
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if  V\  is  of  the  form 

Vt  =  UXm. 

By  substituting  (2.34)  and  (2.35)  into  (2.30)— (2.33),  it  is  easy  to  show  that 


a(X)  = 


*  2  U 

1  +  m 


X'+m,  b(X)  = 


U(\  +  m) 


and  the  boundary  value  problem  then  reduces  to  the  following  simpler  form: 

/"'  +  //"  +  0  {l  -  fn  +  (»  -  1)(1  -  fir)}  =  0, 
g"  +  fg'  =  0, 

/(0)  =  /'(o)  =  0(0)  =  0,  • 

f'(oo)  =  0(00)  =  1, 

where  the  constant  0  and  the  cooling  coefficient  n  in  (2.37)  are  given  by 


0  = 


2m 

1  +  m’ 


n  = 


Tr 


(2.35) 

(2.36) 


(2.37) 

(2.38) 

(2.39) 

(2.40) 

(2.41) 


On  substituting  (2.34a)  into  (2.26a),  we  obtain  u  =  ])  and  insertion  of  this  expres¬ 

sion  and  (2.34b)  into  (2.23)  then  yields  the  following  expression  for  T : 

f  =  r,/«  +  (1  +  i(7  -  1)M!]  {1  -  r  +  (»  -  1 )( I  -  »)}  •  (2.42) 


Our  previous  experience  with  the  zero-pressure  gradient  problem  shows  that  the  stability 
properties  of  a  hypersonic  boundary  layer  depend  very  much  on  the  properties  of  the  basic 
state  temperature  distribution.  When  there  is  no  pressure  gradient,  it  is  known  from  Hall 
and  Fu  (1989)  that  the  boundary  layer  splits  into  two  sublayers  in  the  large  Mach  number 
limit:  a  wall  layer  where  T  =  0(M2)  and  a  logarithmically  thin  temperature  adjustment 
layer  where  T  =  0(1),  and  it  is  the  latter  sublayer  which  supports  the  most  rapidly  growing 
mode.  In  order  to  determine  the  precise  boundary  layer  structure  for  the  present  non-zero 
pressure  gradient  problem,  we  first  need  to  know  the  large  7  behaviour  of  f  and  g.  It  can 
be  shown  from  (2.37)  and  (2.38)  that  as  7  — »  00, 


fin)  =  1  + 


1(1 

\2( 


-  l)c2 


Cl 


2(7  —  ijj)  (7  —  u>y+2(3 


je-^— )a  +  ..., 


g  =  1 - — — e"  +  . .  ■ ,  (2.43) 

7  —  LO 

where  u>,  c\  and  c2  are  constants.  On  substituting  (2.43)  into  (2.42),  we  find  that  as  7  — ♦  00, 

?  =  T>  +  +  •  •  •  •  <2'44) 
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where 


N  =  c  1(7-1). 


In  order  to  determine  the  location  and  thickness  of  the  temperature  adjustment  layer  at 
which  T  =  0(1),  it  is  convenient  to  let 


Tx* 


(ri-uy+W  N 


Solving  this  equation  for  (r/  —  cu)  gives 


T]  —  u>  =  Mi  — 


y  +  (1  -f  2(3)  log  Mi 

Mi 


(2.45) 


where  _ 

Mi  =  \f2\ogAP. 

Equation  (2.45)  shows  that  the  temperature  adjustment  layer  is  at  rj  =  0(M  1)  and  is  of 
0(1 /Mi)  thickness.  In  the  adjustment  layer,  T  is  given  by 


f  =  Ti(  1+e5),  (2.46) 

so  that  T  — ►  Ti  as  y  — »  —00  and  matches  with  the  wall  layer  solution  (2.44)  as  y  — ♦  +00. 

In  summary,  the  basic  state  under  consideration  corresponds  to  Vi  —  Ui/\/T\  =  U Xm 
and  has  solution  given  by 

a  =  i=-^{ptf/1C(x)a\X)f(v)+UJ\1)(Y,+  ^ir2Y)y  (2.47) 

w=W,g(l),  r  =  r1/'2  +  (l  +  i(7-l)MJ){l-/'2  +  /l1(l-s)},  (2.48) 

where  the  new  variables  X,  Y,  Y  and  ij  are  related  to  (x,y)  by 

X  =  f-,M2P,{x)jTMC(x)Jx,  Y  =  yft\(7)Y,  Y  =  £pdy,  r,  =  (2.49) 


and  it  can  be  shown  that  the  function  Yx  appearing  in  (2.47b)  is  given  by 


V  p'  (  TI,  VI n-  1  [  <>nx'Y>dV 

Y-  =  KTL  '  —g^dY- 


(2.50) 


In  section  4,  we  shall  need  to  know  the  asymptotic  expressions  for  these  basic  state  quantities 
and  their  various  derivatives  in  the  temperature  adjustment  layer  in  order  to  determine  the 
relative  orders  of  the  perturbation  quantities.  These  expressions  are  listed  in  the  Appendix 
at  the  end  of  this  paper. 
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3  The  linear  Gortler  instability  problem 

In  order  to  determine  the  evolutionary  behaviour  of  Gortler  vortices  in  the  hypersonic  bound¬ 
ary  layer  described  above,  we  now  superimpose  a  vortex  structure  on  the  basic  state.  Thus 
the  total  flow  is  written  as 


u  =  u  - f  Ue,az,  v  =  v  +  W“2,  to  =  to  -f  We'02, 


p  —  Pi  +  —  Peiaz,  T  =  T  +  0eiaz,  p  =  p  +  peiaz.  (3.1) 

a 

Here  U,  V,  etc.  are  functions  of  X  and  Y  only.  On  substituting  (3.1)  into  the  governing 
equations  (2.9)  — (2.13)  and  linearizing,  we  obtain  the  following  perturbation  equations: 

jvHT)  -  jfL(p)  ~  f(“*  +  »»)  -  f(L(&)  +  ia™e) 


du  dV 


dP  ,  WdP 


+  ^  +  ^  +  r  {/^  +  I/f  =0, 


dx  dy 


dx  dy 


p{L(U )  +  iawU  +  Uux  +  Fuj,}  —  j;L(u)0 
d  t_9U\  .  d 


,  _  2  dU  p  d2V  .  dU  .  pdW 

(pv,  +  pGnu)U  +  -ps- a--— 

2  4  d  (  dV 

+ia  •  -/XyW  +  P  {LiY)  +  iawV)  +  pvyV  +  a2pV  -  -  —  f  p— 

+  3^  -  ^Tx  -  ,C ^  +  -  3 n)Ty 

~  {  +  jGsu2]  -  jO.) v  +  +  (/•«»)«)  9  -  o, 

jxdf/  _  .  /xdK  4_  2 

xu/x*!/  +  xa  •  -  pwxU  +  iapyV  +  ta  ■  -  -pa  W 

d  (  <9W\  „  .  „  .  2 . .  _  .  _  dO 

d’dy  I  /x-^—  1  -  pu>y  K  -  xa/3  -  xa  •  -;x(txr  + 

+  [j;L(w)  -I-  {pwy)y]0  -  p[L(W)  +  iawW)  =  0, 

( PT,  -  (7  -  ')M2^-]G  -  2(7  -  1)M!  •  +  pf.V 

+^(i(«)  +  i<iu>«]  +  o2  ^ 
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(3.6) 


-  +  (7  -  l)A/2/iK)2  +  9  =  0, 


Here  y,  =  /x(T),  fi  =  dy,/dT,  and  the  operator  L(  )  is  defined  by 


(3.7) 


...  _  d  _  d 

i()  =  “&+tV 


In  these  equations,  variables  x  and  y  have  not  yet  been  transformed  into  X  and  Y  using 
(2.49).  This  will  be  done  at  a  1  t‘er  stage. 

So  far  we  have  not  specified  the  size  of  the  wavenumber  a  relative  to  the  Mach  number. 
We  know  from  the  analysis  in  §2  that  the  hypersonic  boundary  layer  splits  into  two  sublayers 
in  the  large  Mach  number  limits:  a  wall  layer  of  0(M\)  thickness  where  T  =  0(M2)  and 
a  temperature  adjustment  layer  of  0(1)  thickness  sitting  at  the  edge  of  the  boundary  layer 
where  T  =  0(1).  We  note  from  (2.49)  that 


dy  = 


b(X) 


Tdrj. 


Hence  in  terms  of  the  physical  variable  y,  the  wall  layer  and  the  temperature  adjustment  layer 
are  of  thickness  of  0(M2Mi)  and  0(1 /Mt),  respectively.  Thus  Gortler  vortices  trapped  in 
the  wall  layer  have  large  wavelength  with  a  =  0(1/{M2M\));  whilst  Gortler  vortices  trapped 
in  the  temperature  adjustment  layer  have  shorter  wavelength  and  have 

a  =  0(±-).  (3.8) 

It  can  be  shown  that  the  growth  rates  for  these  two  modes  are  respectively  of  order  one  and 
order  \/ M2M\.  Thus  the  mode  trapped  in  the  temperature  adjustment  layer  is  the  most 
dangerous  mode  and  will  be  studied  in  the  rest  of  this  paper. 


4  The  inviscid  mode  trapped  in  the  temperature  ad¬ 
justment  layer 

As  in  the  case  when  there  is  no  pressure  gradient  (see  Hall  and  Fu  (1989)),  an  important 
term  in  the  perturbation  equations  is  that  involving  L(v)  in  the  y-momentum  equation.  It 
is  easy  to  show  with  the  aid  of  the  basic  state  properties  listed  in  the  Appendix  that 

L(v)  =  +  (4.1) 
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where  the  expression  for  B(X )  is  given  in  the  Appendix.  This  term  requires  the  Gortler 
number  to  be  of  order  M 2  (the  wall  curvature  should  first  of  all  overcome  the  effective 
negative  curvature  of  the  basic  state  if  Gortler  vortices  were  to  exist  at  all).  Hence  we  scale 
G  by  writing 

G  =  M2G.  (4.2) 


With  the  order  of  G  fixed,  we  can  now  use  the  basic  flow  properties  to  deduce  the  relative 
orders  of  the  perturbation  quantities.  It  is  found  that 

A  .  O(v^),  V  =  0(W)  =  0(^0),  P  =  0(^-6),  (4.3) 

and  that  U  is  of  such  an  order  that  it  does  not  enter  the  leading  order  analysis.  The  large 
growth  rate  implied  by  (4.3a)  is  driven  by  the  term  L(v)  +  \kGu2  in  the  y-momentum 
equation.  This  makes  the  viscous  terms  in  the  perturbation  equations  unimportant  and  so 
the  resulting  eigenvalue  problem  is  essentially  inviscid. 

With  the  use  of  (4.1)  and  (4.2),  we  have 

L{v)  +  ^kGu2  =  Qk(x)G  -  B(x)j  M2. 


Later  we  will  show  that  if  G  is  chosen  such  that  the  coefficient  of  M 2  in  the  above  equation 
vanishes  and  if  the  crossflow  is  zero,  the  0(y/ M2 Mi)  growth  rate  will  become  zero  and 
i«:(.r)(j'  =  B(x )  then  defines  the  leading  order  neutral  position  (this  is  not  true,  however, 
when  the  crossflow  is  non-zero,  as  we  shall  show  later).  In  the  neighbourhood  of  this  leading 
order  neutral  position,  Gortler  vortices  have  a  smaller  growth  rate.  Such  near- neutral  modes 
are  dominated  by  viscous  effects  (see  Hall  and  Fu  (1989)). 

We  now  turn  to  the  determination  of  the  order  of  crossflow  amplitude.  It  can  be  seen 
from  (3.2)— (3.7)  that  this  is  achieved  by  investigating  the  terms  involving  the  operator 

d  d 

u—  +  v—  +  taw. 
ox  oy 

With  the  use  of  the  asymptotic  expressions  listed  in  the  Appendix,  this  operator  can  be 
written  in  the  temperature  adjustment  layer  as 

+  (4.4) 

to  leading  order.  It  may  appear  that  crossflow  effects  will  be  felt  when  d/dX  ~  aW^,  but  it 
will  be  shown  later  that  the  leading  order  constant  term  in  w  does  not  affect  the  leading  order 
growth  rate,  see  Hall  (1985)  for  a  similar  result  for  the  incompressible  problem.  Therefore, 
the  appropriate  size  of  the  crossflow  is  fixed  by  letting 

a  . 
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and  in  view  of  (3.8)  and  (4.3a),  this  gives 


Wx 


M3 


We  thus  introduce  an  0(1)  parameter  6  by 


Wx 


N  M 3 

ciT‘  M?+v 


(4.5) 


to  characterize  the  scaled  crossflow  amplitude.  With  the  order  of  the  crossflow  fixed,  we  now 
look  for  the  following  form  of  WKB  solutions  for  (3.2)-(3.7): 


IM 2 


k  =  ^-(kfl(X,!/)  +  -)£  + 


w  =  \[jr,{W°{x'i)+'")E+'"' 

«  =  (»„(  A •,»)  +  •••)£  +  •••.  (4.6) 

p  =  w,mx'  »)+■")£+•”. 

where  E,  defined  by 

£ (v®^. f  mrwdx  -  f  S  •  <4-7> 

takes  care  of  the  fast  variation  of  the  perturbation  quantities  in  the  streamwise  direction. 
Here  A(X)  is  in  general  complex,  the  real  part  of  it  characterizes  the  growth  rate  of  the 
disturbance.  In  each  of  the  expansions  in  (4.6),  the  first  ellipse  represents  the  higher  order 
terms  which  tends  to  zero  as  M  — >  oo,  whilst  the  second  ellipse  represents  the  meanflow 
correction  and  higher  harmonics.  We  note  that  the  larger  constant  term  iaWx  in  (4.4)  has 
been  absorbed  by  inserting  the  second  term  in  the  expression  for  E.  Thus,  as  remarked 
below  (4.4),  the  term  iaW\  in  w  only  contribute  a  sinusoidal  coefficient  to  the  disturbance 
without  affecting  the  growth  rate. 

On  substituting  (4.6)  into  the  linear  perturbation  equations  (3.2)— (3.7)  and  make  use  of 
the  basic  state  properties  listed  in  the  Appendix,  we  obtain,  by  keeping  only  those  leading 
order  terms, 

~TW  +  ikw' °  =  0,  (4‘8) 

C  \  BP 

(A  -  ikSe^Vo  -  =  0,  (4.9) 

T  px  ay 
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(4.10) 


-ft  +  Sfc'Vi  +  S(A  -  iWe«)lV„  =  0, 
/>1  J  11  1 

-  2Wo  +  (A  -  iMe^o  =  0. 


where 


i, _  6  a  *  _  Pi  V7T  / 1  ^ 

Piv^TA/i’  6  V2° 


(i(jK  -  B(X))  . 


(4.11) 


(4.12) 


It  can  be  seen  that  the  function  G  defined  above  represents  a  balance  between  the  positive 
wall  curvature  and  the  negative  curvature  due  to  the  basic  state. 

By  eliminating  Wo,0o  and  Pq  from  (4.8)— (4.11)  in  favour  of  Vo  and  introducing  a  new 
independent  variable  (  by 

C  =  e*,  '  (4.13) 


we  obtain 


where 


-20!Vo  ,  C(i-0«ft 


C’aTT  +  \~Tr~^F  +  "(Oft  =  0, 


'  acJ  l+c  sc 

mn  >W  C(l-C)  .2 
H{<)  -  ' ITT  + 


\  (A  -  ii 


5-0+02 


(4.14) 


(4.15) 


A 

It  is  convenient  to  eliminate  the  first  order  derivative  term  from  (4.14)  by  introducing  V 
through 


Equation  (4.14)  then  reduces  to 


C2—  + 

<>  ac2 


p  -  r£v°- 


.  C(i-C)  „,n 
+  (l  +  C)2+  ( 


Jv  =  o, 


(4.16) 


where  the  ‘hat’  over  V  has  been  dropped  to  simplify  notation.  Under  the  transformation 
C  =  e*,  — oo  <  y  <  oo  becomes  0  <  C  <  oo  and  therefore  equation  (4.16)  is  to  be  solved 
subject  to  the  asymptotic  boundary  conditions  V  — *  0  as  (  — *  0,  oo.  These  conditions 
should  be  made  more  precise  for  numerical  calculation  and  it  can  be  deduced  from  (4.16) 
that  the  correct  decay  behaviour  is 


5+‘ { 1  ~ StTT11  +  T  +  “  2)1  +  0(<2)}  “ 


C  ke  k<",  as  (  — >  oo. 


(4.17) 


The  eigenvalue  problem  (4.16)  and  (4.17)  determines  the  relationship  between  A  (the  real 
part  of  which  is  the  growth  rate),  the  scaled  crossflow  amplitude  S  and  the  scaled  wavenumber 

A 

k.  In  addition,  G,  defined  by  (4.12),  can  also  be  taken  a s  a  parameter  for  fixed  X.  We  note 
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that  for  each  set  of  fixed  values  of  8,  k  and  G ,  if  (A,  V)  is  an  eigenvalue  set  of  (4.16)  then  so 
is  (  —  A*,  V*),  where  an  asterisk  denotes  complex  conjugate.  Therefore  to  each  unstable  mode 
there  corresponds  a  stable  mode.  Also,  for  each  set  of  fixed  values  of  k  and  G,  if  (6,  A,  K) 
is  an  eigenvalue  set  of  (4.16)  then  so  is  (—5,  —A,  V).  Hence  we  can  restrict  our  attention  to 
positive  crossflow  parameters  and  will  do  so  for  the  remainder  of  this  paper. 

Before  we  solve  this  eigenvalue  problem  for  the  non-zero  crossflow  cases,  let  us  first  note 
that  when  there  is  no  crossflow,  8  =  0  and  then  (4.14)  can  be  written  as 


g  f  C  dVo] 

dC  1(1  +  C)2  J 


k 2 


CU  +  C)2 


j-(i  +  C)!  +  jjcJ  Vo  - 


(4.18) 


Equation  (4.18)  subject  to  Vo  vanishing  at  £  =  0,  oo  is  a  standard  Sturm-Liouville 
problem  which  can  be  taken  to  be  an  eigenvalue  problem  for  G/ A2  for  each  fixed  k.  Its 
solution  is  of  the  form 


G 

A* 


=  F(Jfe),  i.e.  A2  = 


G 


F(kY 


(4.19) 


where  F(k)  is  a  positive  function  of  the  scaled  wavenumber  k  and  can  easily  be  determined 
by  solving  (4.18)  numerically  using  the  fourth  order  Runge-Kutta  method.  This  solution 
will  be  discussed  further  in  the  following  section. 

It  follows  from  (4.19)  that 


A2  >  0  if  G  >  0, 


A2  <  0  if  G  <  0. 


If  we  denote  the  location  where  G  —  0  by  Xn,  then  A  is  real  on  one  side  of  Xn  and  is  purely 
imaginary  on  the  other  side.  This  implies  that  as  X  =  Xn  is  crossed,  the  behaviour  of  the 
disturbance  changes  from  being  oscillatory  to  being  exponentially  growing.  Thus  we  may 
define  X  =  Xn  as  the  neutral  position  and  in  view  of  the  definitions  (4.2)  and  (4.12b),  the 
corresponding  neutral  Gortler  number  G  then  expands  as 


G  = 


2  B(Xn) 
K(Xn) 


M2  +  higher  order  terms, 


(4.20) 


where  the  higher  order  terms  can  be  obtained  by  following  the  procedure  given  in  Hall  and 
Fu  (1989).  However,  our  purpose  in  the  present  paper  is  not  to  determine  the  expression 
for  the  neutral  Gortler  number,  but  instead  we  wish  to  investigate  how  the  existence  of 
a  crossflow  would  affect  the  growth  rate  of  the  inviscid  mode  described  by  the  eigenvalue 
problem  (4.16).  To  be  more  specific,  one  of  our  aims  in  this  paper  is  to  determine  how 
the  growth  rate—  wavenumber  relation  shown  in  (4.19)  for  zero  crossflow  is  modified  as  the 
crossflow  is  gradually  increased. 
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5  Numerical  results  for  non- neutral  modes 


The  eigenvalue  problem  (4.16)  and  (4.17)  has  been  solved  by  using  the  fourth-order  Runge- 
Kutta  method  to  determine  the  dependence  of  the  growth  rate  (the  real  part  of  A)  on  the 
wavenumber  k  with  8  fixed  or  on  the  crossflow  amplitude  8  with  k  fixed.  The  boundary 
C  =  0  is  chosen  to  be  some  small  number  and  (,  =  oo  is  replaced  by  a  large  number  which  is 
continuously  adjusted  for  varying  k  and  8  values  (since  as  we  will  see  later  the  eigenfunctions 
become  more  and  more  concentrated  around  the  singular  point  as  we  approach  neutral 
stability).  At  each  end,  V  is  chosen  arbitrarily  and  dV/dC  is  then  calculated  from  (4.17). 
We  shoot  from  the  two  ends  and  match  the  solutions  at  an  appropriate  middle  point  where 
we  impose  the  condition  that  (\/V)dV/d(  has  the  same  valpe  from  either  side. 

For  increasingly  small  wavenumbers,  the  eigenfunctions  become  more  and  more  spread 
into  the  whole  boundary  layer  (as  we  would  expect  since  as  k  — ►  0  we  are  approaching  the 
wall  layer  mode)  and  in  particular  the  eigenfunctions  vary  rapidly  near  the  wall.  In  this  case, 
it  is  found  convenient  to  introduce  a  scaled  variable  £  through  (  =  e*  for  the  calculations  in 
the  interval  to  the  left  of  the  matching  point. 

In  order  to  be  able  to  continue  the  calculations  even  when  the  mode  becomes  neutrally 
stable  (i.e.  when  A  becomes  purely  imaginary),  we  also  introduce  an  indentation  in  our 
integration  path.  In  order  to  see  whether  the  indented  path  should  be  above  or  below  the 
real  axis,  we  first  note  from  (4.15)  that  the  complex  singular  point  of  the  governing  equation 
is  at  (*  given  by 

C  =  -  A,i),  (5-1) 

where  Ar  +iA,  =  A.  For  growing  Gortler  vortices  Ar  >  0  and  so  this  singular  point  lies  below 
the  real  axis.  Thus  if  the  neutrally  stable  mode  is  taken  to  be  approached  from  the  unstable 
mode,  then  the  integration  path  should  be  indented  from  above.  In  our  calculations,  the 
indentation  is  chosen  to  be  a  triangle  and  checks  were  duly  made  to  make  sure  that  it  does 
not  make  any  difference  whether  the  indentation  is  on  the  left  or  right  hand  side  of  the 
matching  point  and  how  large  the  indentation  is. 

We  now  present  our  numerical  results  for  the  non-neutral  modes.  All  of  the  following 
results  correspond  to  G  =  2  and  are  for  the  first  (i.e.  the  fundamental)  mode.  Other  choices 
of  positive  G  values  can  be  shown  to  yield  qualitatively  similar  results.  The  cases  G  <  0 
will  be  discussed  in  the  final  section.  We  note  that  although  G  <  0  is  a  possibility,  cases 
corresponding  to  G  >  0  are  more  important  since  when  8  =  0,  positive  values  of  G  give  rise 
to  unstable  Gortler  vortices  over  the  whole  wavenumber  regime,  whilst  if  G  <  0,  no  unstable 
Gortler  vortices  can  exist. 

We  shall  first  investigate  the  dependence  of  the  growth  rate  on  the  wavenumber  with 
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8  fixed.  In  Fig.  2(a)  we  have  shown  such  growth  rate  curves  for  a  number  of  8  values, 
including  the  zero-cross  flow  case  8  =  0.  We  see  that  for  8  =  0,  the  growth  rate  tends  to 
zero  as  k  — *  0  and  tends  to  a  constant  as  k  — ►  oo.  The  former  behaviour  is  required  as  the 
present  mode  should  match  with  the  wall  layer  mode  as  k  — ►  0  and  the  wall  layer  mode  has 
smaller  growth  rate.  We  now  show  that  the  latter  behaviour  is  governed  by  a  WKB  solution 
of  the  governing  equation  (4.18). 

When  6  =  0,  Gortler  vortices  are  governed  by  equation  (4.18).  As  k  — ♦  oo,  we  look  for 
WKB  solutions  of  the  form 

Vo  ~  exp (k  JC  E(()d(). 

Equation  (4.18)  requires  that 

C!£2(0  =  fc-(i  +  0!- 

It  then  follows  that  in  general  there  will  be  two  first  order  turning  points  where  E(()  =  0. 
But  the  first  mode  which  we  are  looking  for  must  correspond  to  these  two  first  order  turning 
points  coalescing  to  form  a  single  second  order  turning  point.  This  occurs  when 

(£"2)2"4  =  0’ 

which  yields  the  only  acceptable  result 


The  corresponding  second  order  turning  point  is  then  at  (  =  1.  In  order  to  determine 
the  higher  order  correction  terms  to  the  growth  rate,  we  have  to  carry  out  the  standard 
asymptotic  analysis  for  second  order  turning  point  problems.  Such  an  analysis  shows  that 
in  the  zero  crossflow  case  A  expands  as 

a  =  2  "  +  W3'2')-  (52) 

The  dotted  line  shown  in  Fig.  2(a)  corresponds  to  this  asymptotic  result.  We  note  that  the 
asymptotic  analysis  leading  to  (5.2)  is  very  similar  to  that  detailed  in  Fu,  Hall  and  Blackaby 
(1992)  for  Sutherland’s  law  fluids  and  that  the  growth  rate  curves  for  the  two  viscosity  laws 
have  the  same  structure  but  they  are  fundamentally  different  from  the  growth  rate  curve  for 
incompressible  Gortler  vortices.  We  also  note  that  the  growth  rate  curve  for  8  =  0  given  in 
Fig.  2(a)  is  also  consistent  with  those  growth  rate  curves,  obtained  by  Dando  and  Seddougui 
(1992),  which  emerged  as  the  Mach  number  is  gradually  increased. 
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As  8  is  increased  from  zero  gradually,  we  see  from  Fig.  2(a)  that  the  growth  rate  is 
decreased  at  each  fixed  wavenumber  and  it  no  longer  tends  to  a  constant  as  k  — ►  oo  (this 
behaviour  is  also  seen  in  Dando’s  (1992)  0(1)  Mach  number  calculations).  Instead,  we  see  for 
8  =  0. 1  that  there  exists  a  critical  wavenumber,  say,  at  which  the  growth  rate  vanishes. 
Such  critical  wavenumbers  also  exist  for  other  8  values  (as  can  be  seen  from  the  exact  solution 
to  be  discussed  in  the  following  section)  although  we  are  unable  to  show  their  existence  in 
Fig.  2(a)  due  to  numerical  difficulty.  For  k  >  no  solution  can  be  found  (since  if  a  stable 
mode  could  be  found,  then  an  unstable  mode  could  also  exist).  Thus  at  each  8  value,  Gortler 
vortices  are  unstable  in  the  wavenumber  regime  (0,^),  stable  in  the  wavenumber  regime 
( kCT,oo )  and  are  neutrally  stable  at  k  =  kCT. 

Fig.  2(b)  shows  seven  more  neutral  curves  for  the  larger*  8  values  indicated.  We  see  that 
whilst  for  small  8  values  crossflow  effects  are  stabilizing  over  the  whole  wavenumber  regime, 
this  is  no  longer  the  case  for  large  enough  8  values.  Fig.  2(b)  shows  that  crossflow  effects  are 
actually  destabilizing  over  the  small  wavenumber  regime  whilst  remaining  stabilizing  over 
the  large  wavenumber  regime.  We  note  that  as  8  — *  oo,  the  critical  wavenumber  at  which 
Ar  =  0  tends  to  a  definite  limit.  We  shall  show  in  the  next  section  that  this  limiting  value  is 
given  by  1  /4.  Hence  we  can  deduce  that  for  large  8  values,  unstable  Gortler  vortices  can  only 
exist  in  the  wavenumber  regime  (0,  *)  and  that  for  any  Gortler  vortex  with  a  wavenumber 
lying  in  the  region  (j,oo),  there  always  exists  a  8  value  at  which  the  vortex  is  completely 
stabilized. 

For  Gortler  vortices  with  wavenumber  lying  in  the  region  (0,  A),  numerical  calculations 
show  that  A  =  0(8),  as  8  — ♦  oo.  Therefore  in  the  large  crossflow  limit,  we  can  look  for  the 
following  asymptotic  solutions  for  (4.16)  and  (4.17): 

A  =  Xq8  +  Aj  +  A2^  1  +  •  •  •  5 


v/(Ar,o  =  K,(Ar,o  +  r,K1(x,c)  +  ---. 


(5.3) 


By  substituting  (5.3)  into  (4.16)  and  (4.17),  it  is  easy  to  show  that  A0  and  are  determined 
by  the  eigenvalue  problem 


|Wo  ,  [  1 
dc  1 4 


+ 


C(1-C) 


+ 


ik 


C(1-C)  ,2 


(1  +<)2  Ao  -i*c  1  +  C 


-*'(1+0 


!|  Vo  =  o, 


Vo-C? 


j 


'-TTTk[x  +  l-2e]  +  0(C)^  “  c^0’ 


'}■ 


(5.4) 


Vo  — *  C  *e  fc<.  as  (  -♦  oo. 

We  see  that  this  reduced  eigenvalue  problem  is  now  independent  of  G  and  hence  the  Gortler 
vortices  in  the  wavenumber  regime  (0,  ^)  can  appropriately  be  identified  as  crossflow  vortices 
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since  in  the  large  crossflow  limit,  their  existence  is  dictated  by  the  crossflow  instead  of  by  the 
wall  curvature  effects.  The  eigenvalue  problem  (5.4)  can  be  solved  numerically  to  determine 
the  dependence  of  Ao  on  the  wavenumber  k,  but  we  shall  not  present  the  corresponding 
results  here  since  these  results  merely  confirm  the  large  crossflow  behaviour  of  those  growth 
rate  curves  shown  in  Fig.  2(b). 

In  Figures  3(a,  b),  we  have  shown  the  dependence  of  the  imaginary  part  of  A  on  the 
wavenumber  for  the  same  8  values  as  those  used  in  Figures  2(a,  b).  We  can  see  that  At-  >  0 
for  Ar  >  0,  which  implies  that  when  Ar  =  0,  the  critical  point  defined  by  (5.1)  lies  in  the 
interval  of  integration  and  thus  the  neutral  mode  (to  be  studied  in  the  following  section)  has 
a  critical  layer  structure. 

In  order  to  see  the  dependence  of  the  growth  rate  on  the  crossflow  for  fixed  wavenumbers, 
we  have  shown  in  Figures  4(a,b)  the  growth  rate  curves  in  the  (Ar,6)  plane  for  a  selection  of 
k  values.  We  can  immediately  draw  the  following  conclusions.  For  each  fixed  k  value  greater 
than  1/4,  there  exists  an  unique  critical  8  value,  8„  say,  such  that  Ar  >  0  for  0  <  8  <  8„, 
Ar  =  0  at  8  =  8cr,  and  there  is  no  solution  for  8  >  8 For  k  <  1/4,  we  always  have  Ar  >  0, 
which  confirms  the  findings  shown  in  Fig.  2(b)  that  Gortler  vortices  with  wavenumber  lying 
in  (0,  ^)  are  always  unstable.  For  k  =  1/4,  neutral  stability  is  attained  at  8 „  =  oo.  Figures 
5(a,b)  show  the  dependence  of  the  imaginary  part  of  A  on  the  crossflow  for  the  same  k 
values  used  in  Figures  4(a,b).  Here  again  we  have  A,  >  0.  In  Figures  6(a,b),  we  have  shown 
the  evolution  of  the  real  and  imaginary  parts  of  V  for  fixed  values  8  =  0.2,  G  =  2  as  k  is 
increased.  We  can  see  that  a  critical  layer  structure  is  gradually  emerging  as  A:  is  increased 
towards  the  neutral  value  kcr. 


6  The  neutrally  stable  mode 


Numerical  results  presented  in  the  previous  section  suggest  that  to  each  8  value,  there  cor¬ 
responds  a  critical  k  value  at  which  neutral  stability  occurs  and  for  each  k  value  greater 
than  1/4,  there  again  exists  a  critical  8  value  at  which  neutral  stability  is  achieved.  In  this 
section,  we  shall  determine  the  relation  between  k  and  8  when  neutral  stability  occurs. 

When  Gortler  vortices  are  neutrally  stable,  we  have  Ar  =  0  and  so  we  may  write 

A  =  iA,  (6.1) 


where  A  is  real  and  from  our  previous  numerical  results  we  know  it  is  also  positive.  In  the 
present  neutral  case,  the  governing  equation  (4.16)  may  be  written  as 


r^v  fi  c(i-o  i  cn-0 

^  dC  l4  +  (i+02  c-c  i+C 


cc 


*2(C  -  Ce)2 


-*2(i  +  o2  K  =  o, 


) 


(6.2) 
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where 


<•  -  s  >  «•  (63) 

This  equation  has  a  regular  singular  point  at  C  =  Cc  and  in  the  neighbourhood  of  this  singular 
point  the  two  Frobenius  solutions  can  easily  be  shown  to  be  given  by 

K  - 14  - (c  -  {' +  -  cJ)(c  - « + '  ■ }  ■  (fU> 


where 


l  +  ^  =  x  1  + 


In  general,  the  solution  in  the  neighbourhood  of  the  singular  point  would  be  given  by  AV+  + 
BV. L,  a  linear  combination  of  V+  and  VI,  but  for  the  present  neutral  problem,  we  can  show 
that  either  A  or  B  must  be  zero. 

To  the  above  end,  we  use  Mile’s  (1961)  approach  and  consider  a  function  t(X,  ()  defined 

by 

t(X,0  =  (V'V%  (6.6) 

where  a  subscript  i  signifies  taking  the  imaginary  part  and  a  prime  denotes  partial  differen¬ 
tiation  with  respect  to  (.  It  is  straightforward  to  show  with  the  aid  of  (6.2)  that  dr/<?£  =  0. 
Thus  r  can  only  be  a  function  of  X  only.  But  we  know  r(X,0)  =  r(oo)  =  0  and  so 
r(X,O  =  0. 

However,  if 

V  =  AV+  +  BV_,  (6.7) 

we  then  have 

r  =  (ABTV^Vl  +  A*BVy;){.  (6.8) 

As  we  have  argued  in  the  previous  section,  if  the  neutral  mode  is  taken  to  be  the  limit  of 

the  unstable  modes,  then  the  branch  cut  should  be  in  the  +  plane.  We  therefore  have 


Arg(C 


C  -  C  >  o, 

C  -  Cc  <  o, 


where 


Hence  we  have 


«  -  Cc)*0**0  =  1C  - 

o  f  o,  C  -  Cc  >  o, 

1  l,  C-Cc<0. 


VI  =  V'±e-’>(,±‘,)s. 


(6.10) 


(6.11) 
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On  substituting  (6.11)  into  (6.8),  and  making  use  of  the  result  V+V-  —  V+V!_  =  i/,  we  obtain 

r  =  {^BVeil(1+,,)s}i . 

Hence  by  evaluating  the  above  expression  on  the  left  and  right  hand  sides  of  ( c ,  we  obtain 

(AB*)i  =  0  and  (ABmeir^%  =  0, 

which  imply  that  AB *  =  0  whenever  sinru  ^  0  and  so  either  A  or  B  must  vanish  if  v  is  not 
an  integer  (our  numerical  results  show  this  is  also  true  even  when  v  is  not  an  integer). 

Once  we  have  established  this  result,  we  can  determine  whether  the  solution  is  AV+  or 
BV-  with  the  aid  of  the  numerical  solutions  for  the  neutral  mode  (as  the  limit  of  our  non¬ 
neutral  calculations  discussed  in  the  previous  section).  The  numerical  results  for  G  >  0  show 
that  the  eigenfunctions  are  singular  at  the  critical  point  and  thus  we  deduce  from  (6.4)  that 
the  solution  is  always  given  by  BV 1  (since  v  >  1  when  G  >  0  and  V+  is  regular).  This  can 
also  be  verified  by  calculating  the  phase  jumps  across  the  critical  point  of  the  eigenfunctions; 
we  then  find  that  the  phase  jump  is  always  tt(1  —  v)/2  instead  of  7r(l  +  v)/2.  Here  u  is 
calculated  from  (6.5).  In  the  course  of  calculating  these  phase  jumps  for  a  selection  of  k 
values,  we  noticed  the  possibility  u  =  4A:  and  this  was  then  verified  by  calculations  for  a  wide 
range  of  wavenumbers.  This  result  then  points  to  the  possible  existence  of  an  exact  solution 
to  the  neutral  eigenvalue  problem  (6.2).  After  some  manipulation,  such  an  analytical  solution 
is  found  to  indeed  exist  and  is  given  by 

1 

v  =  Ak,  (c  =  6  =  4A 

16G&  =  (16A:2  -  1)62,  (6.12) 

V'  =  BC"(C-Ce)i-“(C  +  ir'e-*c. 

We  note  that  each  component  in  the  product  on  the  right  hand  side  of  (6.12c)  has  a  distinct 
role:  they  in  turn  represent  the  behaviour  of  the  eigenfunction  at  the  four  regular  singular 
points  C  =  0,  Co  —1*  and  oo. 

The  exact  solution  (6.12)  can  explain  quite  a  few  points  about  our  numerical  results. 
Firstly,  we  can  see  from  (6.12d)  that  k  =  1/4  is  an  important  value.  For  any  fixed  G  >  0, 
neutral  states  can  only  be  obtained  for  k  >  1/4  since  (6.1 2d)  does  not  have  a  solution  for 
8  otherwise.  When  k  =  1/4,  the  neutral  state  is  obtained  at  8  =  oo.  Secondly,  for  a  given 
6  value,  there  exists  a  neutral  k  value  (i.e.  a  k„  at  which  Ar  =  0)  and  for  a  given  k  value 
greater  than  1/4,  there  exists  a  neutral  8  value.  All  these  conclusions  certainly  agree  with 
the  numerical  results  presented  in  Figures  2-5.  Finally,  we  see  from  (6.12b)  that  even  for 
moderate  values  of  k,  (c  can  be  quite  small  and  in  these  cases  the  eigenfunctions  become 
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trapped  near  the  wall,  which  give  rise  to  some  numerical  difficulties.  This  is  why  we  have 
not  extended  the  growth  rate  curves  to  the  neutral  points  in  Fig.  2(a)  although  we  have 
done  so  in  Fig.  2(b). 

After  having  succeeded  in  obtaining  the  above  exact  neutral  mode  solution,  we  also 
considered  the  neutral  mode  eigenvalue  problem  formulated,  but  unsolved,  by  Bassom  and 
Hall  (1991).  As  remarked  by  these  authors,  their  equation  (2.14b)  is  just  a  scaled  version  of 
their  equation  (4.10c)  and  it  governs  the  large  wavenumber  inviscid  neutral  mode.  It  suffices 
to  consider  their  equation  (4.10c),  which  can  be  written  as 

<6i3) 

where  ipc  —  — /§/A.  Using  the  same  procedure  as  that  used  in  obtaining  (6.12),  we  found 
that  (6.13)  also  has  an  exact  solution,  given  by 

A  =  -4= -4,  Vo  =  Aip(tp  —  V>c)~1e-^',  (6.14) 

where  A  is  an  arbitrary  constant. 

The  eigenvalue  problem  associated  with  (6.13)  arises  in  Bassom  and  Hall  (1991)  where 
the  re-emergence  of  the  Gortler  vortex  mode  is  considered  at  high  scaled  crossflow  and 
vortex  wavenumbers.  In  that  paper  it  was  found  that  at  sufficiently  high  crossflow  values 
the  Gortler  mode  is  stabilized  over  a  finite  band  of  nonzero  wavenumbers.  This  stable 
band  of  wavenumbers  increases  with  the  crossflow;  its  left-hand  limit  approaches  that  of  the 
stationary  crossflow  vortex  of  Gregory,  Stuart  and  Walker  (1955),  while  its  right-hand  limit 
is  described  by  the  inviscid  eigenvalue  problem  associated  with  (2.14b)  of  Bassom  and  Hall 
(1991).  The  result  (6.13)  shows  that,  in  the  notation  of  Bassom  and  Hall  (1991),  the  upper 
neutral  wavenumber  is  related  to  the  crossflow  parameter  A  by 

A  =  2v/2 a}/2  ( /r1/2G1/2.  (6.15) 

\  ™  (0) / 

Here  R  and  G  are  the  Reynolds  and  Gortler  numbers  whereas  the  unperturbed  boundary 
layer  is  (u,  R~l^2v,  \w).  The  prediction  (6.15)  is  in  excellent  agreement  with  the  numerical 
results  of  Bassom  and  Hall  (1991)  thus  for  example  the  result  (6.15)  differs  from  the  computed 
value  of  a  by  less  than  5%  when  \R}^2G~X^2  =  4.75. 
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7  Conclusion 


So  far,  we  have  consistently  assumed  *hat  G  >  0.  We  have  shown  that  if  6  =  0,  unstable 
Gortler  vortices  can  exist  in  the  whole  wavenumber  regime.  But  as  8  is  increased  gradually 
from  zero,  this  (0,  oo)  unstable  wavenumber  regime  contracts  to  (0,  A)  where  k  is  related  to 
8  by  (6. 12d)  and  tends  to  1/4  as  8  tends  to  infinity.  Any  Gortler  vortex  with  a  wavenumber 
greater  than  1/4  could  be  completely  stabilized  by  a  large  enough  crossflow  in  the  basic 
state.  We  have  also  shown  that  Gortler  vortices  over  the  wavenumber  regime  (0,  j)  are 
always  unstable  and  for  large  crossflows  they  are  actually  destabilized  by  the  crossflow. 
Since  such  Gortler  vortices  for  large  crossflows  are  no  longer  controlled  by  wall  curvature 
effects,  we  have  identified  them  more  appropriately  as  crossflow  vortices. 

To  complete  our  discussions,  we  now  investigate  what  happens  if  G  <  0.  In  the  zero- 
crossflow  case,  result  (4.19)  shows  that  if  G  <  0,  the  inviscid  mode  has  Ar  =  0  and  so  the 
growth  rate  curve  in  the  (k,  Ar)  plane  is  simply  the  k- axis.  We  now  investigate  how  this 
curve  (i.e.  the  fc-axis)  evolves  as  we  increases  8  gradually.  We  note  that  the  neutral  solution 
(6.12)  is  also  valid  for  G  <  0.  When  G  <  0,  we  deduce  from  (6.12d)  that  for  each  8  value, 
there  exists  an  unique  Ar,  kcr  say,  at  which  Ar  =  0,  and  we  must  have  k „  <  1/4.  This  is 
verified  by  the  numerical  results  shown  in  Figure  7(a,b)  where  we  have  shown  two  typical 
growth  rate  curves  corresponding  to  8  =  5,7  and  to  G  =  —2.  It  can  be  seen  that  for  8  >  0 
and  G  <  0,  the  only  unstable  wavenumber  regime  is  now  (0 ,kcr)  with  kcr  <  1/4,  which  is 
even  narrower  than  that  for  G  >  0.  This  is  a  result  of  the  crossflow  trying  to  destabilize 
and  the  negative  resultant  curvature  trying  to  stabilize  these  (crossflow)  vortices.  It  can  be 
deduced  from  (6.12d),  and  it  is  also  confirmed  by  Fig.  7,  that  k „  tends  to  1/4  from  below  as 
8  — >  oo,  in  contrast  to  the  situation  for  G  >  0  where  kcr  tends  to  1/4  from  above  as  8  — >  oo. 

Figures  8(a.b)  show  that  a  crossflow  vortex  at  a  fixed  8  value  and  with  wavenumber 
lying  in  the  region  (0,  j)  can  be  completely  stabilized  by  decreasing  G,  and  so  although  the 
crossflow  vortices  in  the  wavenumber  regime  (0,  ^)  are  always  unstable  for  G  >  0,  they  can  be 
nevertheless  completely  destroyed  by  a  larg^  enough  resultant  negative  curvature  (represented 
by  G).  Therefore,  crossflows  stabilize  the  unstable  Gortler  vortices  with  wavenumber  lying 
in  the  region  (1/4,  oo)  and  the  negative  resultant  curvature  stabilizes  the  unstable  crossflow 
vortices  with  wavenumber  lying  the  region  (0, 1/4). 

To  summarize,  we  have  shown  in  this  paper  that  in  the  hypersonic  context,  the  crossflow 
which  can  completely  stabilize  the  inviscid  Gortler  vortices  considered  must  necessarily  be 


of  order 


R~1'2 


M 3 

M%+20 


in  terms  of  the  original  physical  variable.  Within  this  order,  we  have  shown  at  precisely  what 
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crossflow  value  a  Gortler  vortex  with  a  given  wavenumber  will  be  completely  stabilized.  It 
is  seen  that  large  wavenumber  Gortler  vortices  can  be  stabilized  more  easily  than  smaller 
wavenumber  Gortler  vortices.  This  may  imply  that  in  practical  situations  small  wavenumber 
Gortler  vortices  could  be  excited  more  easily  than  Gortler  vortices  with  larger  wavenumbers. 

It  is  also  seen  that  crossflow  effects  on  Gortler  vortices  are  fundamentally  different  for 
incompressible  and  hypersonic  flows.  For  incompressible  flows,  Bassom  and  Hall  (1991) 
showed  that  only  at  sufficiently  large  crossflow  values,  were  Gortler  vortices  stabilized  over 
a  finite  band  of  wavenumbers,  whilst  we  have  shown  in  this  paper  that  for  any  non-zero 
crossflow,  hypersonic  Gortler  vortices  are  stabilized  over  a  semi-infinite  wavenumber  regime 
(kCT>o o),  where  kCT  is  determined  by  the  exact  solution  (6.1 2d)  (with  kcr  identified  with  the 
k  there). 

We  have  also  shown  that  in  order  to  quantify  vortices’  response  to  crossflow  effects,  the 
scaled  wavenumber  regime  may  be  neatly  divided  into  two  regions:  (0,  1/4)  and  (1/4,  oo). 
We  have  identified  the  vortices  wavenumber  lying  in  the  first  region  as  crossflow  vortices 
since  in  the  large  crossflow  limit,  their  existence  is  no  longer  supported  by  curvature  effects 
but  instead  by  the  crossflow.  We  have  also  found  the  exact  solution  to  the  neutrally  stable 
mode.  This  exact  solution  should  be  useful  for  studying  the  nonlinear  development  of  Gortler 
vortices  which  is  now  under  consideration. 
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Appendix:  Basic  state  properties  in  the  temperature  adjustment 

layer 
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Figure  1.  A  typical  example  of  a  three-dimensional  boundary  layer:  a  boundary  layer  over 
a  swept  aerofoil. 
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The  scaled  wavenumber  k 


Figure  2.  Dependence  of  the  growth  rate  on  the  scaled  wavenumber  for  G  =  2.  (a)  For 
small  crossflows  S  =  0,0.06,0.1,0.2,0.7, 1.;  the  dotted  line  is  the  asymptotic  result  for  6  =  0 
given  by  (5.2);  (b)  For  larger  crossflows  6  =  1,2,3, ...,  7. 


I_J _ L 

0.00 


1 - 1 - 1 - 1 - 1 — 

0.10  0.20  0.30 

The  scaled  wavenumber  k 


0.40 


Figure  3.  Dependence  of  the  imaginary  part  of  A  on  the  scaled  wavenumber  for  G  =  2.  (a) 
For  small  crossflows  6  =  0,0.06,0.1,0.2,0.7, 1.;  (b)  For  larger  crossflows  6  =  1,2,3,  ...,7. 


Fig.  5(a) 


Fig.  5(b) 


Figure  5.  Dependence  of  the  imaginary  part  of  A  on  the  crossflow  amplitude  for  (7  =  2.  (a) 
For  k  =  0.5, 0.6, 0.7, 1.2;  (b  )for  k=0.1,  0.15,  0.2,  0.25,  0.3. 
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Fig.  6(b)  Eigenfunctions  for  6  =  0.2,  G=1 


Fig.  6(a)  Eigenfunctions  for  6=0.2,  G=  1 
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Figure  6.  Variation  of  the  eigenfunction  for  G  =  2  and  6  =  0.2  as  k  is  increased  towards 
the  neutral  value,  showing  the  emergence  of  a  critical  layer  structure,  (a)  The  real  part  of 
V;  (b)  the  imaginary  part  of  V. 
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The  imaginary  part  of  X  The  real  part  of  X  (growth  rate) 


Fig.  7(a) 


The  scaled  wavenumber  k 


Fig.  7(b) 


Figure  7.  Dependence  of  (a)  the  growth  rate  and  (b)  the  imaginary  part  of  A  on  the  scaled 
wavenumber  for  G  =  —  2  and  6  =  5, 7. 
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Figure  8.  Dependence  of  (a)  the  growth  rate  and  (b)  the  imaginary  part  of  A  on  G  for 
k  =  0.1  and  6  =  5,7. 
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